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We investigate the optical properties of one-dimensional (ID) dimerized Mott insulators using 
the ID dimerized extended Hubbard model. Numerical calculations and a perturbative analysis 
from the decoupled-dimer limit clarify that there are three relevant classes of charge-transfer 
(CT) states generated by photoexcitation: interdimer CT unbound states, interdimer CT exci- 
ton states, and intradimer CT exciton states. This classification is applied to understanding the 
optical properties of an organic molecular material, l,3,5-trithia-2,4,6-triazapentalenyl (TTTA), 
which is known for its photoinduced transition from the dimerized spin-singlet phase to the 
regular paramagnetic phase. We conclude that the lowest photoexcited state of TTTA is the 
interdimer CT exciton state and the second lowest state is the intradimer CT exciton state. 

KEYWORDS: strongly correlated electron systems, spin-Peierls transition, photoinduced phase transition, 
optical conductivity 



1. Introduction 

Currently, there is extensive interest in photoinduced 
phenomena of strongly correlated electron systems. ^"^^ 
These phenomena show drastic changes in magnetic or 
electronic properties triggered by photoirradiation, and 
often accompany subsequent ultrafast relaxation to equi- 
librium states. To reveal the fundamental physics of these 
phenomena is one of important challenges in condensed 
matter physics, and furthermore it would be of great help 
for the application of related materials to new optical 
switch devices. 

Theoretical efforts have also been devoted to this field 
recently. -'^^"^^ From the theoretical point of view, the chief 
ingredient to understand the photoinduced phenomena 
is to clarify the nature of excited states relevant to the 
photoexcitation. Photoirradiation generates optically al- 
lowed excited states (photoexcited states), introducing 
novel characteristics to the system. Some photoexcited 
state may lead the system to an unconventional macro- 
scopic phase that is quite different from any thermody- 
namic one. Thus, to identify the relevant photoexcited 
states is the starting point of the theoretical investiga- 
tion for how the photoinduced phenomena occur. 

In this paper, we present a simple picture to describe 
the photoexcited states of a particular class of strongly 
correlated electron systems, one-dimensional (ID) dimer- 
ized Mott insulators. This picture is derived from a per- 
turbative treatment from the decoupling limit of the 
ID dimerized extended Hubbard model and its numer- 
ical calculations. We demonstrate that there are three 
types of photoexcited charge-transfer (CT) states: (i) in- 
terdimer CT unbound states, (ii) interdimer CT exci- 
ton states, and (iii) intradimer CT exciton states. This 
theoretical understanding is applied to discussion about 
the experimental results of an organic radical crystal, 
l,3,5-trithia-2,4,6-triazapentalenyl (TTTA).^^ This ma- 
terial is a quasi-lD dimerized Mott insulator below the 
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spin-Peierls transition temperature, and shows a charac- 
teristic photoinduced transition from the dimerized spin- 
singlet phase to the regular paramagnetic phase. ^^'^'^ In 
addition, its optical excitations have not completely been 
understood.^"*' In this work, we demonstrate that the 
lowest photoexcited state of TTTA is the interdimer CT 
exciton state and the intradimer CT exciton has a higher 
excitation energy. 

2. Model 

In this work, we use the ID dimerized extended 
Hubbard model for ID dimerized Mott insulators. The 
Hamiltonian is given by 



n 



N-l 



(1) 



where c\ ^ {ci,a) is the creation (annihilation) operator 

of an electron with spin a at site I, and n/^o- = cj ^ci^a- 
The nearest- neighbor transfer integrals alternate a,st{l + 
S),t{l — 6), ■ ■ ■ , and the parameter t, the transfer integral 
of the regular lattice, is used as a unit of energy in this 
paper. The band-filling of electrons is set to be a half. 
The Coulomb repulsion is taken into account up to the 
nearest neighbor: the on-site Coulomb interaction U and 
the nearest-neighbor interaction V. 



3. Analysis from the Decoupling Limit 

In this section, we invoke the analysis from 



the 



decoupled-dimer limit of the model (1)-^^' ^^••^'^ This anal- 
ysis is able to treat U/t exactly, and is proved to be 
appropriate to construct a quite simple picture of the 
photoexcited states. Although several studies have also 
presented different approximations to discuss the opti- 
cal properties of this model, ^^"■^^ these preceding theo- 
ries are justified only in the [//t — > oo limit. In fact, 
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our treatment also has a disadvantage; the vahdity of 
the approximation is not confirmed for small 5. How- 
ever, in this work, we can restrict ourselves to systems 
with rather large dimerization, typically with 6{< 0.4), 
because the organic material TTTA is considered to have 
a quite large value of 5 ~ 0.46.^^ 

In our preceding paper, this analysis for V = has 
provided a clear classification of the photoexcited charge- 
transfer (CT) states: intradimer CT states and inter- 
dimer CT states. In the following, wc discuss how this 
classification is influenced by finite V. All the notations 
used here are defined in ref. 26, with a replacement from 
U toU' = U -V. 

3.1 Model 

Following the formulation in ref. 26, wc split the 
Hamiltonian (1) into the intradimer part Ho and the in- 
terdimer part Hi as follows: 

n = Ho+Hu (2) 

with 

Ho = -t{l +S)J2 J2^4+l,a':i,'^ + 4,aCl+h<^) 
T lee 

CT lee 

+VY,nini+i (3) 

lee 

and 

Hi = -t{l -S)Y, Y,^c}+i,a(^i,'7 + 4ctC;+1,ct) 
CT ;eo 

+Vj2nini+u (4) 

leo 

where X]iee(o) denotes the summation over even (odd) 
/. In what follows, Hq is considered to be the 0-th order 
Hamiltonian and Hi be the perturbation term. Hence 
this treatment is valid for {1 — S),V/t « 1. 

Before discussing photoexcited CT states, we identify 
the ground state of the 0-th order system. The ground 
state l^o) of ^0 is a direct product of the ground states 
of isolated dimers, 

IV'g) = 1(5)0 ®|G)i0-- -0 1(9) ;v,-i, (5) 

where the state |G)„ is the ground state of the n-th dimcr 
(n = 0, 1, • • • ,Nd-l) and Nd{= N/2) is the number of 
dimers. We note that for a strongly correlated system 
([/' >> t), \G)n is equal to the so-called spin-singlet 
dimer (| Ti) - | iT))/V2. 

3.2 Intradimer CT state and interd,imer CT staie 
The intradimer CT state is generated by an intradimer 

charge transfer. By this charge transfer on the dimer n, 
the dimer ground state |G)„ is converted to |0^)„ defined 
by 

The corresponding excitation encrg\ is given by 

Aintra = U' /2 + ^/U^JIT^ , (7) 
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where t' = t{l + S). We note that the state |0^)„ has 
an exciton-like character because the holon and the dou- 
blon^^ are bound in a single dimer. 

The interdimer CT state is generated by an in- 
terdimer charge transfer. By this charge transfer be- 
tween the dimers n and n + 1, the two dimer ground 
states \G)n\G)n+i are converted to |0^)„|E^)„+i or 
\El)n\Ol)n+i. Here, |0^)„ and |£'^)„ are eigenstates of 
the isolated dimer n and defined by 

I'^ct)" — "^('^2nCTC2rs+lCT'^2nCT + '^2n+lCTC2n+lCT<^2nCT)|0)> 
\El)n = ^(4^+4+1^)10), (8) 

where a and a are spin indices and a is opposite to a. 
Its excitation energy is given by 

Ainter = V - 2t' + 2^/U^/I+^. (9) 

In contrast to the intradimer CT state, the interdimer 

CT state has a free-carrier-like character for V/t << 1. 
For V = 0, \0^) and \E^) move freely by the Ist-order 
process of Hi .^^ 

For V/t << 1, the frcc-carricr-like interdimer CT state 
is the lowest photoexcited state and the exciton-like in- 
tradimer CT state has a higher excitation energy.^^'^^'^*^ 
For U' » t', the two excitation energies are given by 

Aintra ~ f/ - V, (10) 

and 

Winter 

U-2t{l + S). (11) 

Thus, the reduction in the energy of the interdimer CT 
state is caused by the kinetic energy gain due to the holon 
in |£^) and the doublon in |0|). 

3.3 Bound-unbound transition of the interdimer CT 
state 

Here we demonstrate that small but finite V causes 
a bound-unbound transition of the interdimer CT state. 
As shown in the previous subsection, the lowest pho- 
toexcited state for V/t << 1 is the interdimer CT un- 
bound state that consists of free "particles" \0^) and 
\E^). Thus, it is reasonable that, for small but finite V/t, 
\0^) and \E^) are bound, resulting in an interdimer CT 
exciton state. To discuss this point, we here use the fol- 
lowing basis set:^^ 

|ni,n2) = ^ [|0?£;J)„,,„, + |Oj4)„„„, 

~I^T*-'i)«l,"2 - \ElO^)ni,n2] ■ (12) 

Here the state \XY)ni^„2, with {X,Y) = {Ol,El) or 
{E^, Og), means that \X) is located on the ni-th dimer, 
\Y) lies on the n2-th {ni < nz) dimcr, and the other 
dimers are in the singlet-dimer state \G). Without H^, 
|ni,n2) are degenerate for arbitrary pairs of (ni < 712). 
Then wc carry out the degenerate perturbation theory 
to treat the first-order effect of 7i^, and obtain 

Hi\ni,n2) =i ^ (|ni + r, 712) -h |ni, 712 -h r)) (13) 

r=-l,l 
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for n2 — ni > 2, and 
Wi|ni,n2) = i{\ni-l,n2) + \ni,n2+l))-V\ni,n2) (14) 
for n2 — ni = 1, where 

i = t{l-S){a + l3f/4: ,a = 2t'/C, , 

P = {^U'^/A + Af^-U'/2)/C, 

C = + 8t'^ - ^7VC^'V4 + 4^'2)l/^(15) 

and 

V = V/4. (16) 

This Hamiltonian is equivalent to that of a two-spinless- 
fermion system defined by 

n = '^[i{ajai+i + aj_^_^ai) - Vnim+i], (17) 

and the "particles" move with the transfer integral t and 
attract with each other via V (see Fig. 1). The attractive 
interaction V is caused by the attractive force between 
the holon in \E^) and the doublon in \0^). The Bethe 
ansatz^^ of eq. (17) tells us that the bound- unbound 
transition occurs a,t V = 2i, resulting in 

Vk/t = 2{l-S){a + (3f. (18) 

Therefore for V < Vb(> Vh) the lowest state is an un- 
bound state (a bound state), and the bound state is the 
interdimer CT exciton state. 



t V t 

^|03>..|E1>-^ 



jG^ ^ ^ jG^ 

Fig. 1. Schematic picture of "particles" in the interdimer CT 
state. 



3.4 Summary of the photoexcited CT states 

Our discussion has demonstrated that there are three 
types of photoexcited CT states: 

(1) The interdimer CT unbound state: when V < Vb, 
\0^) and generated by the interdimer charge- 
transfer, move freely in the whole system. 

(2) The interdimer CT exciton state: when V larger 
than Vfe is introduced, |0^) and \E^) are bound. 

(3) The intradimer CT exciton state: the holon and the 
doublon are strongly bound in a single dimer. 

These three types of CT states have been mentioned in 
a model with longer-range repulsive interactions, associ- 
ated with nonlinear optical responses of conjugated poly- 
mers. Here we demonstrate that these three states are 
important for the systematic understanding of the pho- 
toexcited states in molecular crystals including TTTA, 
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where models with short-range interactions are appro- 
priate.^'' In addition, we show that their wave functions 
and relative energies are governed by V/t. 

Although this classification is basically valid for small 
1 — 5, V/t where the perturbative treatment of Tii is ap- 
propriate, for large V , a simple speculation can be ob- 
tained as below. As V increases, the intradimer CT exci- 
ton state linearly decreases its excitation energy. On the 
contrary, the interdimer CT state is less sensitive to V. 
As a result, a lc;vcl crossing would occur between them 
and then the lowest photoexcited state would change 
from the interdimer CT state to the intradimer CT state. 
In fact, there is a crossover of the lowest photoexcited 
state, from the interdimer CT state to the intradimer 
CT state because two CT states have the mixing term,^*^ 

n{OX+i{G\'Hi\0'')n\E^)n+i = t{l - 5){a + 13). (19) 

This argument is summarized in a schematic energy di- 
agram of Fig. 2. When V <Vb, the lowest photoexcited 
state is the interdimer CT unbound state with \0^) and 
l-E^) while the intradimer CT exciton has a higher exci- 
tation energy. As for the spectral weight, as shown later, 
the intradimer CT exciton is dominant and the inter- 
dimer CT unbound state has a small contribution. As V 
increases, the interdimer CT state becomes the exciton 
state with bound \0^) and \E^), although its spectral 
weight would remain small. As V further increases, the 
intradimer CT exciton lowers its energy level and ap- 
proaches the interdimer CT exciton state, resulting in 
the crossover between them. Then the two states have 
comparable spectral weights because of the mixing. For 
much larger V, the lowest photoexcited state becomes 
the intradimer CT exciton state and has the dominant 
spectral weight. 




V 

Fig. 2. Energy diagram of relevant excited states. The solid (dot- 
ted) lines show the exciton (unbound) states, and the widths of 
the lines indicate the strength of the spectral weight. 



4. Numerical Results 

4-1 Optical conductivity 

In this section, we examine whether our analytical dis- 
cussion is appropriate or not for the realistic parameter 
region. To find photoexcited states of the model (1), we 
calculate the optical conductivity spectrum of the ground 
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state iV'o) given by 

^('^) ^^T. \{^^\J\^o)\^S{uJ + Eo- E,), (20) 



where 



J = [1 + {-l)'5]{c\+,^,cu. - 4,Q+i,,) (21) 

is the current operator, the i-th eigenstate, and Ei 
the corresponding energy. The lattice constant is set to 
be unity for simplicity. We note that the Drude part is 
omitted because the model (1) is treated only in the in- 
sulating phase. In the actual computation, Lorentzians 
with the broadening parameter e are substituted for the 
5 functions in eq. (20), and the Lanczos diagonalization 
method is used to obtain a{uj) and other quantities for 
finite systems {N < 14) with the open boundary. 

Figure 3 shows cr(w) and the excitation energies of the 
two pronounced peaks for 5 = 0.4, U/t = 10, and N = 12. 
Although 8 is rather small in this case, our argument is 
confirmed to be appropriate. For V = Q, there is a dom- 
inant peak at uo/t ^ 11, which corresponds to the in- 
tradimer CT exciton state. As V increases, the excitation 
energy of this peak decreases linearly, and approaches 
that of the lowest photoexcited state ji/'iopt)- It has been 
shown that IV'iopt) for small V is the interdimer CT state, 
and its spectral weight is small. However, as the energy 
of iV'iopt) is approached by that of the intradimer CT ex- 
citon, the spectral weight increases. Around V/t = 2, the 
gap between the two levels is the smallest, and for larger 
V/t the lowest state IV'iopt) has the dominant spectral 
weight, implying the intradimer CT exciton state. 



where Ahd is the projection operator onto the subspace 
that has only one holon-doublon pair, and is the 

number operator of doublon (holon) defined by 



ni^i^ni^l and = (1 - n;^|)(l - n;^|). (23) 

Here we put the doublon at the (A^/2 — l)-th site, which 
is the left side of the central bond connecting two dimers. 
Figure 4 shows Ph(0 for S = 0.4. For V/t — 0, the holon 
spreads over the whole system, refiecting the unbound 
nature of the state. When V/t is introduced, -Ph(6) in- 
creases, and becomes the largest in the whole system at 
V/t = 1.5. Thus, l^/'iopt) is considered to be the inter- 
dimer CT exciton state, where the binding occurs be- 
tween the holon and the doublon that are generated by 
the interdimer CT excitation. As V/t further increases, 
the state IV'iopt) becomes the intradimer CT exciton 
state where the holon and the doublon lie in a single 
dimer. 



V/t=0.0 

v/t=i .5; 

V/t=3.0 




Fig. 4. (Color online) Holon distribution function for S 
0.4, U/t = 10 and N = 12. 




Fig. 3. (Color online) Optical conductivity spectra (solid lines) 
and excitation energies (dotted lines) of photoexcited states for 
5 = 0.4, U/t = 10 and N = 12. The broadening e is set to 0.3t. 



4-2 Holon distribution 

To demonstrate the change of the lowest photoexcited 
state |'0iopt) more clearly, we calculate a quantity associ- 
ated with the nature of the photoexcited state. Here we 
consider the holon distribution function Ph of |'0iopt), 
which is defined by 



Phil) ^ 



_ ('/'lopt|Ahd("-^/2_i'^z')-Plhd|V'l opt/ 



l^'ihdlV'iopt)! 



(22) 



5. Photoexcited States of TTTA 

Let us apply our theoretical conclusion to understand- 
ing the photoexcited states of an organic quasi- ID dimer- 
ized Mott insulator, TTTA.-^^ The optical properties of 
TTTA have been studied experimentally, and it has been 
shown that 17(0;) of TTTA has a double-peak structure. 
However, the origin of the peak structure has not com- 
pletely been understood. "'^'^'^^ 

To address this issue, we first estimate the physical pa- 
rameters of TTTA. As for the dimerization S, 6 = 0.46 
is obtained on the basis of the extended Hiickel calcula- 
tion in ref. 28. The other parameters are determined by 
comparison between the experimental and the numerical 
results of a{uj), and they turn out to be U/t = 10, V/t — 
1.8, and t = 0.235 eV with the Lorentzian width e ~ Q.7t 
(see the dotted line in Fig. 5). The numerical result well 
reproduces the double peak around 1.8eV and 2.1eV. We 
also note that the obtained Coulomb interaction U — 
2.35eV is comparable to that by a first-principles calcula- 
tion {U = 2.9eV),'^* supporting our estimation. For bet- 
ter fitting of the shoulder structure around 2.5eV, slightly 
smaller S = 0.4 and U/t = 10, V/t ^ 1.7, t = 0.235 eV 
with e = 0.7t is also a good parameter set (see the solid 
line in Fig. 5). 
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1 2 CO (eV)3 4 



Fig. 5. (Color online) Comparison between the measured optical 
conductivity of TTTA from ref. 28 (filled circles) and the calcu- 
lated optical conductivity of = 14 systems with the parameter 
sets described in the text (solid and dotted lines). 



Next, we show that the lowest photoexcited state of 
TTTA is the interdimer CT exciton state, and the sec- 
ond peak is caused by the intradimer CT exciton state. 
To iUustrate this conclusion, we calculate the holon dis- 
tribution Ph(0 of the correction vector^^ 

1^) = ^T^T^T^^s— T'^l'^o), (24) 

where A is the normalization factor. This formulation 
facilitates us obtaining photoexcited states with arbi- 
trary CO. Using eq. (24), we obtain Ph{l) of photoexcited 
states corresponding to the two peaks, \uj — 1.8eV) and 
\u! = 2.1eV), and the results for the latter parameter 
set (5 = 0.4, [//t = 10,V/t = 1.7) are shown in Fig. 6. 
The results clearly show that |1.8eV) is an interdimer 
CT exciton state and |2.1eV) is an intradimer CT exci- 
ton state. This assignment is consistent with the larger 
spectral weight of |2.1eV). 




Fig. 6. (Color online) Holon distribution functions in two pho- 
toexcited states |1.8eV) and |2.1eV) for <5 = 0.4, U/t = 10, V/t = 
1.7 and TV = 12 with e = 0.054. 



From the experimental point of view, Fujita et al. 
have made the same assignment of these peaks; |1.8eV) 
is the interdimer CT exciton state and |2.1eV) the in- 
tradimer CT exciton state. This conclusion is derived 
from the dominant spectral intensity of |2.1eV). How- 
ever, the same research group made the opposite assign- 
ment; |1.8eV) is the intradimer CT exciton state and 
|2.1eV) the interdimer CT exciton state. This was de- 
rived from the excitation energies. They have erroneously 



considered that the intradimer CT exciton would have 
the larger binding energy because the holon and the dou- 
blon are confined in a small region. 

Our numerical results support the former conclusion. 
In addition, our analysis from the decoupling limit proves 
that the excitation energy is not always a good measure 
to distinguish the nature of the photoexcited state in this 
system. As shown in eqs. (10) and (11), the excitation 
energy of the interdimer CT state can be lower than that 
of the intradimer CT state for small V/t. It is caused by 
the kinetic energy gain of the interdimer CT state, and 
this is the case for TTTA. For large V/t, the intradimer 
CT exciton state becomes the lowest photoexcited state. 

6. Summary 

We have studied the properties of the charge-transfer 
excitations in the ID dimerized extended Hubbard 
model. Numerical results combined with the analysis 
from the decoupling limit have proved that there are 
three classes of photoexcited states: the interdimer CT 
unbound state, the interdimer CT exciton state, and 
the intradimer CT exciton state. Which state is the 
lowest-energy excitation depends on the nearest-neighbor 
Coulomb interaction V . This classification of the pho- 
toexcited states is applied to the interpretation of the 
optical conductivity spectrum of the organic material 
TTTA. We conclude that the lowest photoexcited state 
(|1.8eV)) of TTTA is the interdimer CT exciton state 
and the higher excited state (|2.1eV)) is the intradimer 
CT exciton state. 
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